Nucleon-Nucleon Coincidence Spectra in the Non— Mesonic Weak Decay of 

A— Hypernuclei and the T n /T p Puzzle 

G. Garbarino, A. Parreno and A. Ramos 
Departament d'Estructura i Constituents de la Materia, Universitat de Barcelona, E-08028 Barcelona, Spain 

(February 8, 2008) 



m 

O 
O 

(N 
C 

a 



The main open problem in the physics of A-hypernuclei is the lack of a sound theoretical in- 
terpretation of the large experimental values for the ratio F n /T p = T(An — > nn)/T(Ap — > np). 
To approach the problem, we have incorporated a AN — > nN one-meson-exchange model in an 
intranuclear cascade code for a finite nucleus calculation of the nucleon-nucleon (angular and en- 
ergy) coincidence distributions in the non-mesonic weak decay of ^He and j^C hypernuclei. The 
two-nucleon induced decay mechanism, Anp — > nnp, has been taken into account by means of 
the polarization propagator formalism in local density approximation. With respect to single nu- 
cleon spectra studies, the treatment of two-nucleon correlations permits a cleaner and more direct 
determination, from data, of r n /r p . Our results agree with the preliminary coincidence data of 
KEK-E462 for ^He. This allows us to conclude that T n /F p for %Re should be close to 0.46, the 
value predicted by our one-meson-exchange model for this hypernucleus. Unfortunately, to disen- 
tangle one- and two-body induced decay channels from experimental non-mesonic decay events, 
three-nucleon coincidences should be measured. 

PACS numbers: 21.80.+a, 13.30.Eg, 13.75.Ev 
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Since many years, a sound theoretical explanation of 
the large experimental values of the ratio, T n /T p , be- 
tween the neutron- and proton-induced non-mesonic 
decay widths, r(An — > nn) and P(Ap — > np), of A- 
hypernuclei is missing [1,2]. The calculations underes- 
timate the central data for all considered hypernuclei, 
although the large experimental error bars do not allow 
one to reach any definite conclusion. Moreover, in the ex- 
periments performed until now it has not been possible to 
distinguish between nucleons produced by the one-body 
induced and the (non-negligible) two-body induced de- 
cay mechanism, ANN —* nNN. Because of its strong 
tensor component, the one-pion-exchange (OPE) model 
with the AI = 1/2 isospin rule supplies very small ratios, 
typically in the interval 0.05 -f- 0.20. On the contrary, the 
OPE description can reproduce the total non-mesonic 
decay rates observed for light and medium hypernuclei. 
Other interaction mechanisms beyond OPE might then 
be responsible for the overestimation of T p and the un- 
derestimation of r„ . Those which have been studied ex- 
tensively in the literature are the following ones: 1) the 
inclusion in the AA — ► nN transition potential of mesons 
heavier than the pion (also including the exchange of cor- 
related or uncorrelated two-pions) [3-7]; 2) the inclusion 
of interaction terms that explicitly violate the AI =1/2 
rule [1,8,9]; 3) the inclusion of the two-body induced 
decay mechanism [10-12] and 4) the description of the 
short range AN — > nN transition in terms of quark de- 
grees of freedom [13,14], which automatically introduces 
AI = 3/2 contributions. 

Recent progress has been made on the subject. 
(1) On the one hand, a few calculations [5-7,13] with 
AN —> nN transition potentials including heavy-meson- 



exchange and/or direct quark contributions obtained ra- 
tios more in agreement with data, without providing, 
nevertheless, a satisfactory explanation of the puzzle. In 
particular, these calculations found a reduction of the 
proton-induced decay width due to the opposite sign 
of the tensor component of A-exchange with respect to 
the one for 7r-exchange. Moreover, the parity violating 
AiV( 3 5i) — > nN( 3 Pi) transition, which contributes to 
both the n— and p-induced processes, is considerably en- 
hanced by A-exchange and direct quark mechanisms and 
tends to increase r„/r p [6,13]. 

(2) On the other hand, an error in the computer program 
employed in Ref. [15] to evaluate the single nucleon en- 
ergy spectra from non-mesonic decay has been detected 
[16]. It consisted in the underestimation, by a factor ten, 
of the nucleon-nucleon collision probabilities. The cor- 
rection of such an error leads to quite different spectral 
shapes and allows one to extract smaller values of T n /T p 
(which is a free parameter in the polarization propaga- 
tor model of Refs. [15,16]) when a comparison with old 
experimental data for ^ 2 C [17] is done. 

In the light of these recent developments and of new 
experiments [18-20], it is important to develop different 
theoretical approaches and strategies for the determina- 
tion of the r„/r p ratio. In this Letter we present a finite 
nucleus calculation of the nucleon-nucleon coincidence 
distributions in the non-mesonic weak decay of \He and 
^ 2 C hypernuclei. The work is motivated by the fact that, 
unlike the single nucleon spectra, correlation observables 
permit a cleaner and more direct extraction of T n /T p 
from data [1] . An experiment performed very recently at 
KEK [18] has actually measured the angular and energy 
correlations that we discuss in this paper. Some prclimi- 
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nary results of the experiment can already be compared 
with our calculations. 

The one-meson-exchange (OME) weak transition po- 
tential we employ to describe the one-nucleon stimulated 
decays contains the exchange of p, K, K* , w and n 
mesons in addition to the pion [6]. The final state in- 
teractions acting between the two primary nucleons are 
taken into account by using a scattering NN wave func- 
tion from the Lippmann-Schwinger (T-matrix) equation 
obtained with the NSC97f potential [21]. The OME de- 
cay rates predicted by this model are the following ones 
[6]: ri = T„ + T p = 0.32, r„/r p = 0.46 for ^He and 

ri = 0.55, r„/r p = 0.34 for \ 2 c. 

The distributions of primary nucleons emitted in two- 
nucleon induced decays and the ratio T 2 /Ti has been 
determined by the polarization propagator method in lo- 
cal density approximation of Ref. [11]: T 2 /Ti = 0.20 for 
^He and T 2 /T 1 = 0.25 for \ 2 C [22]. 

In their way out of the nucleus, the primary nucle- 
ons, due to collisions with other nucleons, continuously 
change energy, direction and charge. As a consequence, 
secondary nucleons are also emitted. We simulate the 
nucleon propagation inside the residual nucleus with the 
Monte Carlo code of Ref. [16]. 

In Fig. 1 we show the kinetic energy correlation of 
np coincidence pairs emitted in the non-mesonic decay 
of %He. The spectra are normalized per non-mesonic 
weak decay. To facilitate a comparison with experi- 
ments, whose kinetic energy threshold for proton (neu- 
tron) detection is typically of about 30 (10) MeV, and 
to avoid a possible non-realistic behaviour of the in- 
tranuclear cascade simulation for low nucleon energies, 
in all the figures of the paper we required T n , T p > 30 
MeV. A narrow peak is predicted close to the Q-value 
(153 MeV) expected for the proton-induced three-body 
process ^He ^ 3 H + n + p: it is mainly originated by 
the back-to-back kinematics (cos9 np < —0.8). A broad 
peak, predominantly due to Ap — ► np or An — ► nn weak 
transitions followed by the emission of secondary (less 
energetic) nucleons, has been found around 140 MeV for 
cos0„ p > —0.8. The kinetic energy correlation for ^He 
nn pairs (Fig. 2) shows essentially the same structure of 
the np distribution just discussed. 

In Fig. 3, which corresponds to the energy correla- 
tion of ^ 2 C np pairs, the narrow peak appearing at 
T n + T p ~ 155 MeV again gets the dominant contribu- 
tion from back-to-back coincidences. The relevance of 
the nucleon final state interactions (FSI) in ^ 2 C relative 
to ^He can be seen from the second, broader peak ap- 
pearing in the region around 110 MeV for ^ 2 C and 140 
MeV for ^He. This peak is in fact more pronounced for 
the heavier hypernucleus. Another consequence of the 
different FSI effects in ^He and ^ 2 C is the different mag- 
nitude of the tail of the back-to-back distribution at low 
energies. 

Figs. 4 and 5 show the opening angle correlations of 



nn, np and pp pairs emitted in the decay of ^He and \ 2 C, 
respectively. Comparing both figures for nn and np coin- 
cidences, one sees that the back-to-back peaks are more 
pronounced for \He (less sensitive to FSI than \ 2 C) than 
for j^ 2 C, while the (almost uniform) tail of these distribu- 
tions (fecded by FSI) is more significant in ^ 2 C than in 
A He. Since at least one proton of any pp coincidence is a 
secondary particle, the pp spectra are quite uniform: ac- 
tually, due to the relevance of the back-to-back kinemat- 
ics in the weak decay, these distributions slowly decrease 
as cos 9 pp increases. Again as a consequence of FSI, the 
number of pp pairs is considerably larger in ^ 2 C than in 
iHe. 

The ratio r„/r p is defined as the ratio between the 
number of primary weak decay nn and np pairs, 
and N™p. However, due to nucleon FSI, one expects the 
inequality: 
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to be valid in a situation, such as the experimental one, 
in which particular intervals of variability of pair opening 
angle, A9\ 2 , and sum energy, A(Ti+T 2 ), are employed in 
the determination of N nn and N np . Actually, as one can 
deduce from Figs. 1-5, not only N nn and N np but also the 
ratio N nn /N np depends on A8 12 and A(7\ + T 2 ). The 
numbers of nucleon pairs N nn , N np and N pp discussed 
up to now are related to the corresponding quantities for 
the one-nucleon (NjP N ) and two-nucleon (NfP N ) induced 
processes [the former (latter) being normalized per one- 
body (two-body) stimulated non-mesonic weak decay] 
via the following equation: 
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Table I shows the dependence of N nn , N np and 
Nnn/N n p on A9 12 and A(7\ + T 2 ) for ^ 2 C. For com- 
parison, the same quantities for the primary weak decay 
nucleons are listed as well. Without any restriction on 
9nn and the nucleon energies, one notes a great increase 
(by about one order of magnitude) of both the nn and np 
numbers when the effect of the FSI is taken into account. 
Again because of FSI, the use of an energy threshold T^ 1 
of 30 MeV supplies N nn and N np values considerably re- 
duced with respect to the ones obtained at T^ h = MeV. 
On the contrary, the ratio between the number of nn and 
np pairs is much less sensitive to FSI effects. 
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FIG. 1. Kinetic energy correlations of np pairs emitted per 
non-mesonic decay of ^He. See text for details. 
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FIG. 4. Opening angle correlations of nn, np and pp pairs 
emitted per non-mesonic decay of a He. See text for details. 
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FIG. 2. Kinetic energy correlations of nn pairs emitted per 
non-mesonic decay of iHe. See text for details. 
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FIG. 3. Kinetic energy correlations of np pairs emitted per 
non-mesonic decay of A 2 C. See text for details. 
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FIG. 5. Opening angle correlations of nn, np and pp pairs 
emitted per non-mesonic decay of A 2 G See text for details. 

In Tabic II (III), the ratio N nn /N np for ^He (\ 2 C) is 
given for different combinations of opening angle inter- 
val and nucleon energy threshold. In parentheses we also 
report the predictions obtained when the two-nucleon 
induced decay channel is neglected. The results of the 
figures and tables presented in this work are in a form 
that permits a direct comparison with KEK-E462 data 
[18], which arc now under analysis. The first prelimi- 
nary results of KEK-E462 are quoted in Table II. Wc 
note that while the central values of the data are more 
in agreement with the calculation that neglects the effect 
of the two-body induced decay mechanism, the complete 
results are compatible with the upper limits of the same 
data. This could indicate an almost negligible effect of 
the Anp — > nnp process and/or a T n /T p ratio slightly 
lower than the one (0.46) predicted by our OME model 
for %He. To clarify this point, on the one hand one has 
to wait for the final results of KEK-E462. On the other 
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hand, our calculations show that three-nucleon coinci- 
dences are required to disentangle the effects of one- and 
two-body stimulated decay channels from observed de- 
cay events. The simplistic picture that the back-to-back 
kinematics is able to select one-nuclcon induced pro- 
cesses is in fact far from being realistic. 

In conclusion, our OME weak interaction model sup- 
plemented by FSI through an intranuclear cascade sim- 
ulation provides two-nucleon coincidence observables 
which reproduce the preliminary KEK-E462 results for 



^He. This allows us to 
should be close to 0.46. 



conclude that r„/r p for 



iHe 



Although further (theoretical 
and experimental) confirmation is needed, in this paper 
we think we have proved how the study of nucleon co- 
incidence observables can offer a promising possibility to 
solve the longstanding puzzle on the r„/r p ratio. 

In a forthcoming (long) paper [23] we shall discuss the 
nucleon correlation observables for the (one- and two- 
nucleon stimulated) non mcsonic decay of ^He and \ 2 C 
in a sistematic way. Single nucleon spectra will be fur- 
ther subject of this work. In addition, one should treat 
the case of ^H, which is of extreme importance in order 
to test the validity of the AI =1/2 isospin rule in the 
AN — ► nN weak transition [1,9,20], another key point for 
the solution of the r„/r p puzzle. 

This work is partly supported by EURIDICE HPRN- 
CT-2002-00311, by the DGICYT BFM2002-01868, by 
the Generalitat de Catalunya SGR2001-64 and by INFN. 
Discussions with H. Outa are acknowledged. 
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TABLE I. Results for N nn , N np and N nn /N np correspond- 
ing to the non-mesonic decay of A 2 C. A null (30 MeV for the 
numbers in parentheses) nucleon energy threshold and two 
different opening angle regions are considered. 





cos Onn < —0.8 


all 9nn 


/V wcl 


0.20 (0.19) 


0.25 (0.24) 


7V wd 

± * Tip 


0.57 (0.56) 


0.75 (0.72) 


n^/n^ = r„/r p 


0.34 (0.34) 


0.34 (0.34) 


N nn 


0.44 (0.11) 


3.15 (0.33) 


N np 


1.05 (0.26) 


8.40 (0.87) 


N n „/N np 


0.42 (0.43) 


0.38 (0.39) 



TABLE II. Predictions of N nn /N„ p for %Re correspond- 
ing to different nucleon thresholds and pair opening an- 
gles. The numbers in parentheses correspond to calculations 
with T2 = in Eq. (2). The (preliminary) data are from 
KEK-E462 [18]. 



(MeV) 


cos 6n n 

< -0.8 < -0.6 < -0.4 all 


30 


0.61 (0.52) 0.61 (0.51) 0.60 (0.50) 0.54 (0.45) 
0.52 ±0.11 0.50 ±0.10 0.51 ±0.10 


50 


0.63 (0.52) 0.61 (0.51) 0.60 (0.51) 0.56 (0.46) 


TABLE III. Same as in Table II for i 2 C. 


T$ (MeV) 


cos Onn 

< -0.8 < -0.6 < -0.4 all 


30 


0.43 (0.37) 0.43 (0.37) 0.43 (0.37) 0.39 (0.35) 


50 


0.41 (0.35) 0.40 (0.35) 0.40 (0.35) 0.38 (0.34) 
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